Feshbach molecule formation in a Bose— Fermi mixture 
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We investigate magneto-association of ultracold fermionic Feshbach molecules in a mixture of 40 K 
and 87 Rb atoms, where we can create as many as 7 x 10 4 40 K 87 Rb molecules with a conversion 
efficiency as high as 45%. In the perturbative regime, we find that the conversion efficiency depends 
linearly on the density overlap of the two gases, with a slope that matches a parameter-free model 
that uses only the atom masses and the known Feshbach resonance parameters. In the saturated 
regime, we find that the maximum number of Feshbach molecules depends on the atoms' phase- 
space density. At higher temperatures, our measurements agree with a phenomenological model 
that successfully describes the formation of bosonic molecules from either Bose or Fermi gases. 
However, for quantum degenerate atom gas mixtures, we measure significantly fewer molecules than 
this model predicts. 



Weakly bound ultracold Feshbach molecules can be 
created in an atomic gas by sweeping a magnetic field 
across a Feshbach scattering resonance. This process, 
called magnetoassociation, has proven to be one of the 
most robust methods of Feshbach molecule formation and 
is effective for both wide and narrow Feshbach resonances 
as well as for gases of fermionic atoms, bosonic atoms, or 
mixtures of two atomic species Magnetoassociation 
was an essential initial step in the recent success in creat- 
ing a gas of ultracold polar molecules, where the Feshbach 
molecules were subsequently transferred to the ground 
ro- vibrational state Q. This polar molecule experiment 
starts from an ultracold Bose-Fermi gas mixture and cre- 
ates fermionic 40 K 87 Rb molecules; however, the efficiency 
of the magnetoassociation is relatively low § and this is 
a limiting factor in the final molecule number. Fermionic 

Feshbach molecules have also recently been created in 
23 Na _40 K and 6 L j_23 Na m i xtures [J^ In au of triese 

experiments the fraction of atoms converted to molecules 
is less than 30%. 

In contrast, molecule conversion efficiencies as high as 
90% have been observed in homonuclear systems where 
this process has been studied in detail @ . The efficiency 
of magnetoassociation has been shown to depend on the 
magnetic-field sweep rate and on the phase-space den- 
sity of the parent atomic gas. In the limit of sufficiently 
slow sweeps, a phenomenological model, which uses a 
pairing criteria based on the separation of two atoms in 
phase space and has a single fit parameter, was shown to 
describe the measured molecule creation efficiencies for 
different initial gas conditions @. 

Motivated by the relevance to ultracold polar molecule 
experiments as well as an interest in the nature of pair- 
ing in Bose-Fermi mixtures, we report here an investi- 
gation of Feshbach molecule association in an ultracold 
gas mixture of 40 K and 87 Rb atoms. Our measurements 
significantly extend the range of atom conditions studied 
for this mixture 0, @~G3 • We characterize the density- 
dependent timescale of magnetoassociation as well as the 
roles of competing inelastic processes. Guided by these 
results, we are able to produce up to 7 x 10 4 Feshbach 



molecules. 

Our experiments begin with an ultracold mixture of 
87 Rb atoms in the \F, mp) = |1, 1) state and 40 K atoms 
in the |9/2, — 9/2) state, where F is the total atomic 
angular momentum and mp is its projection onto the 
magnetic-field axis. The gas is confined in a far-detuned 
optical trap formed by a single beam of 1090 nm light 
focused to a 19 /im waist. Typical radial and axial trap- 
ping frequencies for our measurements are 450 Hz and 6 
Hz, respectively, for Rb The K trap frequencies are 
higher than those of Rb by a factor of 1.4. We perform 
measurements above the phase transition temperature, 
T c , for Bose-Einstein condensation (BEC) of the 87 Rb 
gas, so that the spatial distributions of the two clouds 
are well-matched. The number of atoms in each species 
is varied between 5 x 10 4 and 6.5 x 10 5 , and the temper- 
ature, T, ranges from 250 nK to 770 nK. 

The atomic samples are prepared at a magnetic field 
that is several Gauss above a broad s-wave Feshbach res- 
onance between the 87 Rb and 40 K atoms. Near the Fes- 
hbach resonance, the scattering length is given by a = 
a bg (1 - A/(B - B )), with B = 546.62 G, A = -3.04 
G, and a& g = —187 a , where a is the Bohr radius 0]. 
We associate atoms into molecules by sweeping the mag- 
netic field downwards across the resonance. One pair of 
magnetic coils is used to provide the bulk of the bias field, 
while a much smaller pair produces fast magnetic-field 
changes on the order of 10 G via controlled discharge of a 
30 mF capacitor. With these coils, and employing open- 
loop corrections for transient magnetic fields caused by 
eddy currents @ , we can realize magnetic-field sweeps as 
fast as 300 G/ms. This allows us to explore the magnetic- 
sweep-rate dependence of molecule formation. 

The remainder of this paper is organized as follows. 
Sec. U describes the details of our magnetoassociation pro- 
cess with a particular focus on minimizing the effects of 
inelastic collisions. In Sec. |TT] we show that the sweep- 
rate dependence of 40 K 87 Rb magnetoassociation follows 
a simple parameter-free model adapted from Landau- 
Zener physics. In Sec. IIIII we measure the saturated 
number of molecules that are produced under a wide 
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range of atom conditions. These are compared against 
the stochastic phase-space sampling model described in 
Ref. @. 



I. MAGNETOASSOCIATION IN THE 
PRESENCE OF INELASTIC LOSS 

The number of molecules detected after magnetoas- 
sociation can be significantly reduced by inelastic colli- 
sions. Molecule-atom collisions, for example, can result 
in vibrational relaxation of the molecule, which releases 
enough energy that both the molecule and atom are lost 
from the trap [lfjj. In addition, three-body inelastic col- 
lisions of atoms can lead to significant inelastic loss and 
heating as the magnetic field approaches a Feshbach reso- 
nance. To minimize these effects, we carefully adjust the 
initial atom gas conditions and our magnetoassociation 
procedure. This facilitates a comparison of our data to 
theories that ignore losses and heating. 

The top part of Fig. Q] illustrates a typical magnetoas- 
sociation sequence. First, the field is swept at a rate of 
3.33 G/ms from 549.9 G down to 545.7 G to form Fesh- 
bach molecules. To selectively image only the molecules, 
we sweep the field in 50 /is down to 544.7 G where un- 
paired 40 K atoms are transferred from the |9/2,— 9/2) 
state to |9/2, —7/2} state using RF Adiabatic Rapid Pas- 
sage (ARP). At this magnetic field, the binding energy 
of the molecules is /ix3 MHz Q, where h is Planck's 
constant, and the ARP for atoms does not affect the 
molecules. To image the molecules, we then sweep the 
field in 65 /is back up to 546.0 G, where the molecules can 
be selectively imaged with light tuned to the cycling tran- 
sition of 40 K. We find that the molecule number measured 
with this imaging technique is 85±5% of the number mea- 
sured after dissociating the molecules (by sweeping the 
magnetic field upwards across the Feshbach resonance) 
and imaging the resulting atoms. Therefore, when imag- 
ing the molecules below the Feshbach resonance, we ap- 
ply a multiplicative factor to correct the number for this 
measurement inefficiency. 

To minimize the loss of molecules due to inelastic col- 
lisions, we can release the gas from the optical trap af- 
ter the molecules are formed. The sudden turn-off of 
the optical trap initiates an expansion of the gas that 
rapidly lowers the densities and switches off inelastic col- 
lisions. Figure [1] shows the measured molecule number 
as a function of the time that the optical trap is turned 
off, Please, where t — is defined as the time when 
the magnetic field crosses the Feshbach resonance. For 
^release > we observe a decrease in the molecule num- 
ber with increasing ircicasc- We attribute this loss to in- 
elastic collisions as the molecules spend more time in a 
relatively high density atom-molecule gas mixture in the 
trap. For t rc \ casc < 0, the magnetic field crosses the Fesh- 
bach resonance after the gas starts expanding. For short 
expansion times, where — 0.3 < Release < 0, the measured 
number of molecules is constant within the measurement 
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FIG. 1. Top: A typical magnetic-field sweep for molecule as- 
sociation. Bottom: Measured molecule number vs the time at 
which the optical trap is suddenly turned off to release the gas, 
Release- The error bars in all figures correspond to plus or mi- 
nus the standard error obtained from repeated measurements. 
Here, t — is defined as the time when B crosses the Feshbach 
resonance. For increasing i rc icasc > 0, the molecules spend 
more time in the trap where atom-molecule collisions can 
cause loss. For comparison, the dashed line shows an exponen- 
tial decay with a time constant of 2.7 ms. For i rc i casc < 0, the 
expansion of the atom gas prior to magnetoassociation even- 
tually reduces the efficiency of molecule creation for a fixed 
magnetic-field sweep rate. For this data, the initial atom gas 
has 3.6 x 10 5 40 K atoms and 1.4 x 10 5 87 Rb atoms at T = 380 
nK (T/Tc = 1.7). 



uncertainty. For longer expansion times, corresponding 
to ircicasc < —0.3, the rapidly dropping atom gas density 
reduces the efficiency of magnetoassociation for the 3.33 
G/ms sweep rate. Thus, we find i re iease = gives the 
maximum number of molecules. 

Another issue is that inelastic collisions can cause 
loss and heating of the atoms even before molecules are 
formed. In particular, three-body recombination is an 
inelastic collision process for three atoms that produces 
a molecule plus an atom, both with excess kinetic energy. 
The rate for three-body recombination increases dramat- 
ically near the Feshbach resonance. For our Bose-Fermi 
mixture, the dominant three-body recombination process 
removes two 87 Rb atoms and one 40 K atom from the 
trap To estimate the size of this effect, we look at 
the number of atoms remaining after a double magnetic- 
field sweep. Figure [2] shows the measured fractional loss 
of 87 Rb atoms as a function of the initial gas density. Af- 
ter magnetoassociation, we dissociate the molecules by 
sweeping the magnetic field back up across the Feshbach 
resonance at a rate of 20 G/ms. Here, the gas is released 
from the optical dipole trap at t = 0. Since the adiabatic 
magnetoassociation process is reversible, this sequence 
should conserve the atom number in the absence of in- 
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FIG. 2. Fractional loss of the 87 Rb atoms after a double 
magnetic-field sweep to create and then dissociate 40 K 87 Rb 
Feshbach molecules. The measured loss is plotted as a func- 
tion of the product of the initial peak 40 K density and the 
peak 87 Rb density squared. Here, the initial atom gas has 
3 x 10 5 S7 Rb atoms and 2 x 10 5 40 K atoms with a Rb radial 
(axial) trap frequency of 490 (6) Hz. The atom densities were 
varied by modulating the optical trap intensity to parametri- 
cally heat the cloud to temperatures between 490 nK and 730 
nK. Inset: magnetic double-sweep sequence. The optical trap 
is turned off at t=0 as the magnetic field crosses the Feshbach 
resonance. 



elastic collisions. Characterizing the density by the peak 
87 Rb density squared times the peak 40 K density, we find 
that the atom loss is below 10% as long as this density 
product remains lower than 2x(10 13 cm" 3 ) 3 . The ini- 
tial atom conditions for data in Sections [TT] and Mil were 
chosen to satisfy this condition. 

Although lower density conditions were chosen to char- 
acterize the molecule formation process, we note that our 
largest molecule clouds are produced from higher atom 
densities. For example, we created 7 x 10 4 40 K 87 Rb 
molecules from a mixture whose density product n^riK 
was 3 x (10 13 cm~ 3 ) 3 . The initial gas in this case had 
6.5 x 10 5 40 K atoms and 2.1 x 10 5 87 Rb atoms at 560 nK. 



II. DEPENDENCE OF 
MAGNETO-ASSOCIATION ON THE 
MAGNETIC-FIELD SWEEP RATE 

The efficiency of molecule creation, /, which is defined 
here as the number of molecules divided by the lesser 
of the initial numbers of atoms of 40 K and 87 Rb, 7V<, 
depends on the magnetic-field sweep rate as shown in 
Fig. [3j Only a small number of molecules are formed 
when the field sweeps too quickly across the Feshbach 
resonance. As the sweep becomes slower and therefore 
more adiabatic, molecule production increases and finally 



FIG. 3. Fraction of the minority atom species converted to 
molecules, /, vs the inverse magnetic-sweep rate, 1/l-B], of 
a 7 G linear ramp of the magnetic field across the Feshbach 
resonance. The dashed line shows a fit to Eq. fT]), which gives 
T = 6.5(5) G/ms and f = 14(2)%, or 3.8(5) x 10 4 molecules, 
for this data. The initial gas consisted of 2.6 x 10 5 40 K atoms 
and 3.4 x 10 s 87 Rb atoms at T= 490 nK. 



saturates. The saturated conversion efficiency for Bose- 
Fermi mixtures is typically well below unity, as is the 
case in Fig. [3l In this section, we focus on the fast-sweep 
regime, where a simple extension of a two-body picture 
predicts that the molecule number increases linearly with 
the sweep duration [l], [l2[ ■ The maximum molecule con- 
version efficiency in the saturated regime is the topic of 
Sec. [ITTJ 

In order to extract information about the fast-sweep 
and saturated regimes independently, we fit the molecule 
conversion efficiency as a function of the magnetic-sweep 
rate \B\ to the following formula: 



/ = fo 1 



-r/(/o|B|) 



(1) 



Here, fo is the saturated molecule conversion fraction and 
/ = r/|B| in the fast-sweep limit. 

In the fast-sweep regime, the molecule number is pre- 
dicted to follow a Landau Zcner-likc behavior dm with 
the transition probability into the molecular state given 
by P — 1 — er 2lx&LZ w For two particles in a 

box of volume V, the Landau-Zener parameter 5lz de- 
pends on the Feshbach resonance parameters and is given 
by El 



SLZ = TT 



1 2irh 



a bg 



A 



D 



(2) 



where h — h/2n and fi is the two-body reduced 
mass. Interestingly, this result applies for thermal atoms 
(fermions or bosons) as well as for atoms in the motional 
ground-state, such as atoms in a Bose-Einstein conden- 
sate [Hj]. Using classical probability theory, Eq. @ can 
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FIG. 4. The initial molecule conversion efficiency divided by 
the magnetic-field sweep rate F vs the density overlap of the 
two atom clouds. The prediction for F using the 40 K- 87 Rb 
Feshbach resonance parameters (Eq. ©) is shown as a solid 
line. Here, the atomic densities are varied by associating 
molecules at different times during expansion from the optical 
trap (up to 1.2 ms), as well as by adjusting the evaporation 
trajectory and optical trap frequencies. For this data, the 
87 Rb radial (axial) trap frequency ranged from 360 to 520 Hz 
(4 to 7 Hz). 



be generalized to multiple particles simply by multiply- 
ing P by the total number of available atom pairs to get 
the number of molecules [l2T |. Using a local density ap- 
proximation, the molecule density at position r is then 
given by 



. . . . , . 2irH 
n mo \(r) = ni(r)n 2 (r)27r 



a bg A 



D 



(3) 



Integrating Eq. ([3]) over the trapped atom gas distribu- 
tion, and dividing by N < , gives 



/ = < 



. 2ttH 
n>)27r 



a bg A 



D 



(4) 



where (n>) = jj— J rti(r)ri2(r)d 3 r describes the density 
overlap of the two atom clouds. 

Plugging the relevant Feshbach resonance parameters 
into Eq. © yields 



r/(n>) = (2.8 ±0.1) x 1CT 12 cm 3 G/ms 



(5) 



This prediction for T, which is shown as a solid line in 
Fig. |4l agrees well with our measurements of T as a func- 
tion of the density overlap of the initial atom gas, (n>). 
From this we conclude that magnetoassociation in the 
perturbative regime of fast magnetic-field sweeps behaves 
as expected in our Bose-Fermi mixture. In the next sec- 
tion, we turn our attention to the saturated regime, which 
is important for experiments where one would like to cre- 
ate larger numbers of Feshbach molecules. 
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FIG. 5. Top: The measured conversion efficiency / vs the ini- 
tial T/Tf of the 40 K atoms. The atom number ratio ranges 
from 1.6 < N K /N Rb < 2 (squares), 2 < N K /N Rb < 4 (cir- 
cles), or 4 < Nk/Nio, < 7 (triangles). The crosses include a 
correction to /, based on a calculation of inelastic collisions 
between 40 K 87 Rb molecules and 87 Rb atoms, shifting / by 
an amount between 0.01 and 0.07. The molecule conversion 
efficiency increases for higher atomic phase-space density or 
lower T/Tf- Bottom: The difference between the SPSS model 
and the data, f ~ fspss, vs the initial T/Tf of the 40 K atoms. 



III. SATURATED MOLECULE NUMBER 

In the saturated regime of slow magnetic-field sweeps, 
molecule creation has been shown to depend on the en- 
tropy or phase-space density of the atom gas (f3,0,[^- A 
corresponding phenomenological model has successfully 
described molecule creation efficiency for homonuclear 
molecules in 40 K and 85 Rb gases as well as in heteronu- 
clear mixtures 85 Rb- 87 Rb [3, 40 K- 87 Rb @, and 6 Li- 40 K 
17]. While previous experimental work with 40 K- 87 Rb 
trapped gas mixtures focused on RF association [1, 0], 
we consider here magnetoassociation. 

Figure [5] shows the measured molecule creation effi- 
ciency / as a function of T/Tp of the 40 K gas, where Tp 
is the Fermi temperature. Here, we varied the magnetic- 
field sweep rate, using the findings of Sec. [TT] to ensure 
that the sweeps were sufficiently slow that the molecule 
number was in the saturated regime. For this data, the 
number of 40 K atoms is larger than the number of 87 Rb 
atoms. Rb radial (axial) trap frequencies range from 
360 to 550 Hz (4 to 6 Hz), and initial atom tempera- 
tures are between 250 nK and 770 nK. As was observed 
for homonuclear molecules @ , we find that the molecule 
conversion efficiency increases for higher initial atomic 
phase-space densities. 
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We compare our results with the Stochastic Phase- 
Space Sampling (SPSS) model that has been used previ- 
ously to treat homonuclear molecule formation 0] ■ This 
semi-classical model is based on the idea that two atoms 
are able to pair if they lie within a volume of relative 
phase-space set by an adjustable parameter 7: 

A* \v re i \ \r re i \ < -^7 (6) 

where v rc \ is the relative speed of the two atoms and r re \ 
is their separation. Another property of the model is that 
atoms are only allowed to pair once (i.e., paired atoms 
are removed from further consideration). To determine 
the parameter 7 in our system, we fit the molecule cre- 
ation efficiency data from Fig. [5] to the SPSS model. We 
first use a Monte-Carlo calculation to generate atom dis- 
tributions matching the initial conditions for each of our 
measurements. These atoms are then paired according 
to the criterion in Eq. ([B]) to determine the molecule cre- 
ation efficiency. A fit to the data with initial T/Tp > 0.55 
returns a value of 7bf=0.38(2). The error bar is domi- 
nated by a 10% systematic uncertainty in the measured 
conversion fraction. Our measured jbf for 40 K- 87 Rb is 
consistent with the value 7^ = 0.38(4) found in 40 K and 
just below the 73 = 0.44(3) measured in 85 Rb [6|. 

In the bottom panel of Fig. [5j we show that the 
SPSS model with 7bf=0.38 significantly overestimates 
the conversion efficiency for atom clouds with low T/Tp. 
We note that previous measurements for homonuclear 
Feshbach molecules in 40 K agreed with SPSS predictions 
for T/T F as low as 0.05 The failure of the SPSS 
model to describe our measurements in the quantum- 
degenerate regime may be due to heating by inelastic 
collisions, or may suggest that this model is inadequate 
for describing the production of fermionic molecules from 
quantum-degenerate atom gas mixtures. 

The maximum molecule conversion efficiency we ob- 
serve in Fig. [5] is 45%, which occurs at our lowest T/Tp 
for the initial 40 K gas. Starting from even lower T/Tp, 
or higher phase space density, would seem then to be 
an obvious way to increase the molecule conversion ef- 
ficiency. However, going to higher initial phase space 
density in the Bose-Fermi mixture would require creat- 
ing a BEC of 87 Rb or having an even larger ratio of the 
number of fermionic atoms to bosonic atoms. In cither 
case, this would result in a reduction in the spatial over- 
lap of the Bose gas and the Fermi gas in the trap, as 
well as in momentum space, which will presumably de- 
crease the molecule conversion efficiency. This effect has 
been observed below the BEC transition temperature T c 
for heteronuclear mixtures [H, [l6[ . In addition, the high 
density of the BEC can result in a rapid loss of molecules 
due to inelastic atom-molecule collisions. 

Another challenge for heteronuclear molecules, such 
as fermionic Feshbach molecules, is that differences in 
the masses, polarizabilities, and/or magnetic moments 
of the initial atoms can cause the two atom clouds to 
have different equilibrium positions in the trap. To mea- 
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FIG. 6. Molecule conversion efficiency, /, as a function of the 
separation between the 40 K and 87 Rb clouds in units of the 
width, o-Rb, of the S7 Rb cloud, which is the smaller of the 
two atom clouds. The solid curve shows an empirical fit to 
a Gaussian with an RMS width of 1.8(1) o-Rb- For this data, 
molecules are created from 2.8 x 10 5 87 Rb atoms and 2.1 x 10 5 
40 K atoms at T = 490 nK and T/T c =1.6. 

sure the significance of this effect on the molecule con- 
version efficiency, we intentionally separated the 40 K and 
87 Rb clouds along the weak axial direction of the op- 
tical trap by applying a magnetic-field gradient before 
molecule formation. The results of this measurement are 
shown in Fig. [SJ where we plot the saturated molecule 
conversion efficiency against the spatial offset of the two 
atom clouds. We find that the conversion efficiency drops 
when the spatial separation is significant compared to 
the root-mean-squared (RMS) size of the 87 Rb gas in the 
trap, o-Rb- For the data in Fig. O the spatial offset of the 
atom clouds in the axial direction is less than 1.2 CRb 
and the calculated offset in the vertical direction due to 
differential gravitational sag is less than 0.3 ORb- These 
relative displacements of the atom clouds were accounted 
for in the SPSS modeling. 

IV. CONCLUSIONS 

We have studied the efficiency of converting atoms 
to fermionic Feshbach molecules in a regime where in- 
elastic collision effects are minimized. We find that in 
the regime of fast magnetic-field sweeps, the molecule 
creation efficiency is proportional to the inverse sweep 
speed with a coefficient that depends only on the Fesh- 
bach resonance parameters and the density overlap of the 
two atom clouds. For sufficiently slow magnetic sweeps, 
where the molecule number saturates, we find that the 
maximum conversion efficiency increases with increasing 
phase-space density. Our measurements with T > T c for 
87 Rb and T/Tp > 0.5 for 40 K agree with a phenomeno- 
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logical model having a single fit parameter. However, atoms than bosonic atoms and for T just above T c . 
significant discrepancies are observed between this model 
and our data for a more quantum-degenerate Fermi gas. 

The highest conversion efficiency, which is about 45% ACKNOWLEDGMENTS 
for our data, is achieved for larger numbers of fermionic 
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the NSF, and NDSEG. 
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